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Quantum Field Theories in Nonextensive Tsallis Statistics 

Hiroaki KOHYAMA*) and Akira NiEGAWA**) 
Graduate School of Science, Osaka City University, Osaka 558-8585, Japan 

Within the framework of TsaUis statistics with q ~ 1, we construct a perturbation theory 
for treating relativistic quantum field systems. We find that there appear initial correlations, 
which do not exist in the Boltzmann-Gibbs statistics. Applying this framework to a quark- 
gluon plasma, we find that the so-called thermal masses of quarks and gluons are smaller 
than in the case of Boltzmann-Gibbs statistics. 

§1. Introduction 

It is expected that a quark-gluon plasma will soon be produced in ultrarelativistic 
heavy-ion collision experiments at the BNL Relativistic Heavy Ion Collider (RHIC). 
The CERN Large Hadron Collider (LHC) will also soon be ready for experiments. 
A few observations regarding quark-gluon plasmas are in order. 

1) According to hot QCD (statistical QCD at high temperature)/^ in the chromo- 
electric sector of gluons in a quark-gluon plasma, a Debye screening mass devel- 
ops, and as a consequence the chromoelectric-gluon exchange interaction is short 
range. By contrast, a Debye-like mass does not appear (at least to one-loop order) 
in the chromomagnetic sector of gluons, and therefore the chromomagnetic-gluon 
exchange interaction is long range. ***^ 

2) The size of a region containing a quark-gluon plasma produced in the collision of 
heavy ions in RHIC experiments is not very large. In fact, the radius of heavy 
ions of mass number ~ 200 is ~ 7.9 x 10"^^ m. At the highest RHIC energy, 
y/SNN = 200 GeV, the Lorentz-contraction factor is ~ 1/100, and therefore the 
longitudinal size of the system just after the collision is Li ~ 7.9 x 10^^'' m, and 
it increases due to the expansion of the system. Contrastingly, the chromoelectric 
Debye mass, given by^-* niD — ^■'^gT (with g the QCD coupling constant and T 
the temperature of the plasma) for three quark flavors, is ~ 2.8 GeV for T = 1 
GeV (where /^n — 0.43), which corresponds to the Debye screening length^) 
Id ~ 6.9 X 10"^'' m. As mentioned above, the chromomagnetic "Debye mass", 
if it exists, is at most of 0{g^T). Thus the "Debye screening length" in the 
chromomagnetic sector is much larger than Iq. 



*' E-mail: kohyamaQsci. osaka-cu.ac.jp 
**■* E-mail: niegawa@sci.osaka-cu.ac.jp 

***-' It is worth mentioning here that, according to hot QED, the magnetic mass does not appear in 
the "transverse-photon" sector to all orders in perturbation theory.'^' Therefore, a magnetic-photon 
exchange interaction in an electron-positron-photon plasma is long range. 

The quark counterpart^' of the Debye screening length is much longer, specifically. Id ~ 
2.1 X 10"^^ m for T = 1 GeV. 
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3) The single-particle transverse momentum distribution of hadrons, which are pro- 
duced by heavy ion collisions, exhibits a power-law tail. 

Boltzmann-Gibbs statistics is valid for the systems that have the following prop- 
erties: (i) compared to the size of the system, both the interactions and the memories 
are of short range, and (ii) the spacetime in which the system evolves is nonfractal. 
Prom the above observations 1) and 2), we find it doubtful that Boltzmann-Gibbs 
statistics can be applied in a strict sense to a quark-gluon plasma, at least in the 
early stages after a collision. It has been argued that the above observation 3) may 
be a sign of the fact that the object (a quark-gluon plasma or a hadronic fire ball) 
produced just after a heavy ion collision does not obey Boltzmann-Gibbs statistics 
(see, e.g., Ref. 3) and references therein). 

For a realistic description of a systems with long-range interactions, long-range 
memories and/or fractal structure, noncxtensive generalization of Boltzmann-Gibbs 
statistics is essential. Such a generalization was proposed by Tsallis^) seventeen 
years ago (so called Tsallis statistics). Since then, a variety of works have appeared 
concerning theoretical aspects of this form as well as its applications to various 
nonextensive systems. The most important problem among the former is to de- 
termine the extent to which Tsallis statistics is unique among other nonextensive 
statistics. A short review of this problem is given in §2.1. Tsallis statistics has 
been successfully applied to a number of nonextensive systems. Examples*^ are 
Levy-type anomalous diffusion,^) Euler turbulence,^^ the specific heat of the hydro- 
gen atom,^) peculiar velocities in galaxies, and self-gravitating systems and related 
matters.'^)' Some theoretical frameworks in Boltzmann-Gibbs statistics have been 
generalized to the case of Tsallis statistics, e.g., linear response theory,^^) the Green 
function method^^^ and path integrals. Some approximation schemes, such as the 
(1 — q) expansion, factorization approximation,^^)' perturbation theory^^^ 
and the semi-classical expansion.^®) 

Tsallis statistics contains a (real) parameter q, which is a measure of nonexten- 
siveness of the system, or more precisely, the range of interactions acting among the 
constituents of the system. For short-range interactions, standard Boltzmann- 
Gibbs statistics is realized, which corresponds to Tsallis statistics with q = 1- 

With the above motivation, in this paper we treat quantum field systems on 
the basis of Tsallis statistics assuming that q differs from 1 by a very small amount. 
More precisely, our region of interest is that satisfying 

\l-q\ <0{l/{VT^f) ^1, (M) 

where T is the "temperature" and V is the volume of the system. **) (For |1 — i^j 
expansions in different contexts, see, e.g., Refs. 14) and 15), together with references 
therein.) 

*■* For a comprehensive list of references, see Ref.5). 

The volume of a heavy ion with mass number ~ 200 is F ~ 2.7 x 10"* (MeV/Rc)~^. (As 
mentioned above, at early stages after the coUision, due to the Lorontz contraction, the volume of 
the produced quark-gluon plasma is much smaller than this V.) For T = 200, 500 and 1000 MeV, 
VT^ ~ 2.1 X 10^,3.3 X 10" and 2.7 x 10^, respectively. 
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In §2, brief introductions to Tsallis statistics and the closed time-path (CTP) 
formalism^)' ^"^'^^^ for treating quantum field systems are given. In §3, we deduce the 
form for the CTP propagators in Tsallis statistics. In §4, we discuss physical impli- 
cations of the results and present a procedure for computing higher-order corrections 
to the propagators. 



§2. Preliminaries 



2.1. Tsallis statistics 

Throughout this paper, we use units in which ks = h = c = 1. Tsallis postu- 
lates'^^ the following form for the generalized entropy (the Tsallis entropy): 

sM = - (2-1) 

q-l 

Here q is a real-number parameter and p is the density operator (Tr p = 1). In the 
limit g ^ 1, Eq. H2-l() reduces to the standard Boltzmann-Gibbs-Shannon entropy 
= — Trplnp. However, contrast to the case of Boltzmann-Gibbs statistics, the 
form for the entropy in the case we consider presently cannot be uniquely deduced. 
The Tsallis entropy H2-l|) meets the requirement of "concavity", the requirement 
which any entropy should satisfy. Since the introduction of the Tsallis entropy, 
dos Santos^^^ has demonstrated the uniqueness of the form appearing in Eq. I\'2-1\\ 
assuming pseudo-additivity, Sa+b = Sa + Sb + — q)SASB (where 'A' and 'S' 
represent the systems in question, and + i?' represents the composite system 
of and '-B'), together with several other conditions. Hotta and Joichi^^^ have 
shown that Eq. H2-l() can be "derived" from less restrictive requirements, namely, 
the composability condition, Sa+b = ^{Sa,Sb) (with 17 being some function), and 
the ansatz S[p] = C + Tr(f){p) (with C a constant and (/) some function of p), together 
with a few other conditions. Modified Tsallis entropies have also been proposed by 
several authors (see, e.g., Refs. 24) - 27)). The connection between Sg and the 
theory of quantum groups has been pointed out and discussed, e.g., in Refs. 26) and 
28). 

The form of p is determined by generalizing the procedure employed in Boltzmann- 
Gibbs statistical mechanics. There had been some disagreement regarding the defi- 
nition of the expectation value of an operator A, but this issue has been settled, and 
it is now known that this expectation value is given by the following:^^^ 

M = — , (2-2) 

which is called the g-expectation value and preserves various desirable properties. 
The form of the density operator p is determined by maximizing Sg [p] with the con- 
straints Trp = 1 and {H) = E, where H is the Hamiltonian. Introducing Lagrange 
multipliers, we easily carry out the maximization and obtain 



1/(1-9) 

l-{l-q)H/T , (2-3) 
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Tr' 



1 - (1 - q)H/T 



1/(1-9) 



(2-4) 



Here, TV' means that Tr is taken over the energy eigenstates with 1 — (1 — q)E/T > 
with E the eigenvalue of H. For 1 — g < 0, this restriction obviously does not apply. 
As mentioned in §1, (7 = 1 corresponds to Boltzmann-Gibbs statistics, p oc e~^^ 
with T = 1/ (3 the temperature. In the following, we simply refer to T in Eq. (|2-4)) 
as the temperature. (For a thorough discussion of the temperature of nonextensive 
systems, see, e.g., Ref. 30) and references therein.) 

Before moving on, for convinience, we rewrite the formula (|2-2|) for the q- 
expectation value as 



{A)=T,'p'A, 



P 



1 - (1 - q)H/T 



9/(1-9) 



1 - eH/T 



{e = l-q) 



1-ePH 



1/^- 



/3^(l-e)/r. 



(2-5) 

(2-6) 
(2-7) 



1-e' 

2.2. Closed time-path formalism 

For treating quantum field systems, we employ the closed time-path (CTP) 
formalism.^)' ^'^^'^^^ In the single-time representation of the CTP formalism, every 
field becomes two fields: (p ~^ {4'1j4'2)- Here, {i = 1,2) is called the type-i field, 
and ^1 is called the physical field. Then, the n-point Green function consists of 2'^ 
components. At the very end of calculation 0i and (j)2 are set equal. 

In the following, we employ the complex scalar field theory governed by the 
Lagrangian density C = —(j)\d'^ + m^)(f) + Cmt with Cmt = —{X/4:){(f>^(f>)'^. Gener- 
alization to other field theories is straightforword. We restrict our consideration to 
the case in which the density operator is electrically neutral. 

2.2.1. Single-time representation 

The 2n-point Green function, which consists of 2^*^ components, is defined by 



G 



n---«n;ii---jri (^1 



,yv 



Tr 



\l=l m=l 

where the operators (f)^^^ and 0^^^^ are the Heisenberg field operators, and p is 
the density operator. Note that, for Tsallis statistics, we have p = p' [Eq. H2-5() ]. 
Here, Tc is the "ordering" operator with the following properties: i) move the type-2 



(2- 



fields to the left of the type-1 fields, ii) rearrange the operators 



*,(-f^) 



according to 



a time-ordering (Tc — > T), and iii) rearrange the operators (p^ according to an 
anti-time-ordering (Tc T). 

Here we summarize the Feynman rules for computing G perturbatively. The 
rules are the same as in the vacuum theory, except that the (bare) propagators and 
vertices take the following forms. 
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1. Propagators (bare two-point functions) 
la) Two-point propagators: 



iAu{x 


-y) 


iA22{x 


-y) 


iAi2{x 


-y) 


iA2i{x 


-y) 



4>i = 4>2 = ' 



(2-9) 



where cj) and cj)^ are the interaction-picture fields and {■■■)= Tr ■■■ po with 
po the "bare" density operator. 

lb) "2n-point propagators" (2 < n): The 2n-point propagators consist of 2^" 
components, all of which are the same: 



Cni^li'"' i^n'iyii ' ' ' lyn) 1 

= {■ 4>(.xi) ■ ■ ■ (i)(xn)(l)Hyi) ■ ■ ■ 4>\yn) ■■) 



,(2<n) (2-10) 



where : • • • : represents the operation of taking take a normal product and 
'c' stands for the contribution from the connected diagrams. Cn is called the 
initial correlation. 

2. Vertices 

The vertex factor for a vertex at which type-1 fields meet is the same as in the 
vacuum theory, i.e. —i\ while the vertex factor at which type-2 fields meet is 
i\. 

2.2.2. Physical representation 

Let us introduce 4>c ^'^id 4>A (and their hermitian conjugates) through the rela- 
tions ^ 

<t>c = -j{<t>l + (t>2) , <t>A = (j)!- (t>2- 

Using these relations, the single-time representation outlined above can be trans- 
formed into the representation written in terms of ^t^A^ (fc and (f)]^, which is called 
the physical representation. 

The 2n-point Green function is defined by*^ 

Gc...cA-A;c-cA-Aixi, - ■ ■ , y„) 

nil n m2 n 



z(-)"Tr 



j=mi+l k=l l=m2+i 



=1 



P , 

(2-11) 



The propagators and vertices in the Feynman rules are as follows. 



*^ The normalization of G here is different from that of its counterpart in Ref. 21). 
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1. Propagators 

la) Two-point propagators: 

iAaa = , 

iA^c = ^(^11 - ^2l) = i'^A , 
iAcA = - A12) = iAr , 



(2-12) 



where ^_r(a) is the retarded (advanced) propagator, which, in momentum 
space, reads 



A 



R(A) 



p2 _ jjj2 _|_ ipQQ+ 



(2-13) 



lb) Initial correlations: Among the 2^" components of the 2n-point propagators 
Cn (2 < n), only the (cc • • • c; cc • • • c)-components are nonzero; 



{xi, • • • • • • ,y„) = C„(xi, • • • • • • ,y„) ,(2-14) 

c--c;c--c 

where C„ is as in Eq. (|2-10|) . 
2. Vertices 

The vertex factors for the vertices 4>l4>l4>c4>A-, 4>c(t>^A'^c^c-, 4>\4>\<pA(kA and (?i'^</>^0c'/'zi 
are —iX/A and — iA/4, respectively. All other vertices vanish. 

2.2.3. Forms of the propagators 

From this point, we restrict our attention to the case po = po{Ho), with Hq the 
free Hamiltonian. Then, the system under consideration is spacetime-translation in- 
variant, and therefore can go to momentum space. To obtain the forms of the prop- 
agators, we first construct single-particle wave functions by adopting a fixed volume 
{V) quantization with discrete momenta, p = (^27r/V^^^) n, where n = (ni,n2,n3), 
with ni, 71-2, ns integers. The complex scalar fields (/'(x) are decomposed by using the 
plane-wave basis constructed in this manner, 

1 



5 i2E^V)y^ 



[ 



(2-15) 
(2-16) 

(2-17) 



The form of Aij{P) is obtained by substituting Eq. (|2T5|) and its hermitian 
conjugate for (p and (j)^ , respectively, into Eq. (|2-9j) . Then, taking the large- volume 
limit, y — > cxD, we obtain (after taking the Fourier transform) 

A,{p) = Af;\p) + A^iP) , 



where Ep = y^p^ -|- is the single-particle energy. Then, Hq becomes 

Ho = Y.Ep(alap + blb^ 
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4?(p) = -(4?(p) 

A^{P) = -2TnN{po)5{P^ - m") 



1 



P"^ -w? + ^0^ 



(2-18) 



where P^ = Pq — and N{pq) is the number-density function, 



Nipo 



hm 



6'(po)TraJ,ap/5o + 6'(-po)Tr6j,6p/5o 
= TriV(po) . 

Similarly, we obtain the form for C„ [Eq. (|2-l()j) ] as follows: 



(2-19) 



C 



n 



d'P; 



n I 



27r6{P^)e{p^)\Y([lN{p,o) n N{-p,o))c 

/ 1=0 i=l j=l+l 

^-iPi-{xi^-yj^) _ _ _ g-^Pri^ii-yji) _|_ perms. I 



L(2vr) 

E E [{ 



iPn-{x., -y, ) 

o n — L ■^n — L 



+ perms. 



(2-20) 



Here, P-{x-y) = po(a;o-yo)-p-(x-y) and X]ii<...<i, (Eji<...<ji) is the summation 
over all possible choices of • • • ii{ji ■ ■ ■ ji) from among the values 1, 2, • • • n, subject 
to the conditins ii < i2 < ••• < ii (ji < J2 < ••• < ji)- i'l 
Ui < j'2 < " ' < j'n~i) obtained from 1, 2, • • • , n by removing ii < i2 < • • • < ii 
(ji < J2 < ■ ■ ■ < ji)- The first 'perms.' indicates that all permutations among 
(ii)j2)''" )Jn) are taken and the second 'perms.' indicates that all permutations 
among (jj, jg, • • • are taken. 

For Tsallis statistics [cf. Eq. (|2-5() and ()2-6l) ]. we have 

lA 



PO 



1 - epHo 



Tr' 



1 - e(3Ho 



l/e ■ 



(2-21) 



This is invariant under charge conjugation, and therefore we have N{pq) = A^(|po|)- 
Equation (|2-2()|) is simplified as 



C 



n 



d^p 



2n6{Pl)\{llN{M)c 

/ 1=1 



(2-22) 
yn) are taken. 



1 {2nY 

X |^e"*-Pi-(^i-2'i) . . . n-v-n-yn; _^ perms. 

where 'perms.' indicates that all permutations among (yi,y2> ■ ■ ■ 
2.2.4. Gibbs ensemble 

A standard Gibbs ensemble with temperature T (= and vanishing chemical 
potential is described by po = e~^^°/Tre~^^°. The initial correlation C„ (2 < n), 
given in Eq. (|2-1U|) . vanishes. 
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From Eq. (|2-17j) with Eq. (|2-16|) . in the limit ^ ^ oo, we obtain 



Tre 



Mm n ( 1 



exp 



-2 X y 



a^(IpoI) 



IPO 



(2^ 

hm ^— 5- — 7 = 

(e/3^^P_l)-i ePoPy 

= iVBE(^p) . 



In ( 1 - e-'^^p 



(2-23) 



,. Tr6Upe-^-f^« 
hm ^—^ — ; 



1 



(2-24) 



Here Vo is the pressure of the free complex scalar field system and A'^be is the familiar 
Bose distribution function. The factor of 2 in Eq. (|2-23|) corresponds to the number 
of degrees of freedom of the complex scalar field. A straightforward manipulation of 
Eq. (IT^ yields, for m/3 << 1, 



Cbe 



g2xCBEW(3/3^) 
7T^ ( . 15 



(2-25) 
(2-26) 



More generally, for a system that consists of single or several kinds of bosons and/or 
fermions, we have 



Tre 



exp 



(2-27) 



where i and j label the kinds of bosons and fermions, respectively. The quantity 
n^j [n^lj) is the number of degrees of freedom of the ith kind of boson (jth kind of 
fermion) and 



•-^BE — 


-1 

30 


(- 


— ( 


r^j) _ 


£! 

30 


a- 


15 
8^ 



(2-28) 



It is worth mentioning that when we replace Tr • • • in Eq. H2-23() with Tr' • • • (for 
e > 0) [cf. Eq. (|MT|l ]. terms of 0(e-^/7e^) and 0{{(3m/efe-'^/^) appear in Eqs. 
(fT^ and ((2I2HI)- Such terms can safely be ignored for e < 1 [cf. Eq. (HHJ]. 



§3. Computation of the propagators 



In this section, we compute the propagators on the basis of Tsallis statistics. 
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3.1. Preliminaries 

~4 



(Tg, appearing in Eq. (|2-21|) . is expanded as follows: 



(Tq = exp 



^ In ( 1 - ipH, 



(3-1) 



e-^^» exp 



-f3Ho 



n 

1=1 



■E 

1=1 

T- 



i+i 



n=0 



l + l 

l + l 



(3-2) 



The region of interest is that defined by Eq. In the following, we compute 

field propagators up to and including 0(e) terms, with e < 0(1/ iVT^Y). For this 
purpose, for (Tq, we should employ the form 



l=Q 



E 

1=0 



21 f2 



2/+4 



+ 



3 



\ 21+2, 



mo] --[PHo) + — [pHo] 



21 



/3 



/3 



18 



/3 



/3 



-/3Ho 



(3-3) 



and for (Tre~^^« ) the form 



/3 



N 



QN 



-/3Ho 



(-)^ 



13 



N 



,2CbeV/(3/33) 



'3N 



{2CBEVf + 2N{N - l){2CBEVy'~'p 



N~l o3 



^ 2N{N -1){N -2){3N -A){2CbeV)''-^ , 



,2CbeV/(3/3») 



(3-4) 



Here, the terms represented by "• • • " yield at most 0(e^/^) contribution to the prop- 
agators. 

In the reminder of the paper, for simplicity, we restrict ourselves to the massless 
case, m = 0, and therefore we have Cbe = 7r^/30 and Ep = |p| = p. Generalization 
to the massive case is straightforward. 

Equation 1)3 •3() and the remark made at the end of §2 guarantee that Tt'ctq ~ 
Tro-g. Then, using Eqs. (EH, and (jHD, we get 
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l + 4{2/(/-l)+/} ^-j 



-1(1 - 1){6(/ - 2)(l - 3) + 23a - 2) + 12}- 
3 y 



-{2l{l - 1) + 7/ + 3}e + -^e - fe + 



2 , , , /lOO o 131 

1 + -y2p(5y - 3) + ey i—f - —y + 24 ) + 



„2CbeV/(3/32) 

(3-5) 



where 



(f = 1//3) 



(3-6) 



3.2. Two-point propagator Aij = 1,2) 

The propagator Aij{P) is given by Eq. H2-18p . with A^(|po|) = N{p) 
the number-density function N{pq). Then, from Eqs. (|2-23|) and H2-24() . we have 



p = \po\ 



for 



hm TraUr,e"^^o = lim Tr5l6r,e"^^° 



_1_^2CbeV/{3P^) = iVBE(p)e2^BEV^/(3/3^) . (3.7) 



Using Eqs. (j^^ with m = and (EHj), we obtain 



lim Tr alcana a 



i iVe2^BEW(3/33) 



1 



- X 1/2 

2Wp{l + N)[-j 



+ 



- {81(1 - l){2l - l)/3p(l + N) + l{2l - 1)CPp?{1 + N){1 + 2N)\ % 
2 1 J y 



-(2/ + 3)/3p(l + iV)ey + 



-2/?p(l + iV) ( fj + 4/?p(l + iV)y(3 - 2y)e 

2 ~ 



1/2 



+ A^(2Ar^ + 3Af + l)(2y - l)e - -/3p(l + iV)y(3 - 2y)e + 



g2CBEV/(3/33) 
(3-8) 
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where = l/{e^P — 1). From Eqs. H3-5|) and H3-8() . we finally obtain for the number- 
density function, 

V^oo IrCJn 



iV 



1 - V2ye/5p(l + iV) 



+£<! 6/3p(l + iV)y + ■^M!(2iv2 + m + l){2y - 1) } + 



Nbe 



2ye/?p(l + iVBE) 



e\{l + 6y)(3p{l + NBE) + 



+0(6^/2) 



(1 + Nbe) (1 + 2iVBE) 



(3-9) 



where use has been made of Eq. (|2-7p . In Eq. (|3-9p . we have used A^'be = l/(e^P — 1) 
and 



■3^2 



(3-10) 



Generalization to other fields is straightforward. For a system of bosons and/or 
fermions, we have, for the number-density function, 



2y-l 



1 - ^/2^ePp{l ± iVBE/FD) + + 6?/)/?p(l ± A^be/fd) 



+- 



{(3pY{l ± iVBE/FD)(l ± 2iVBE/FD) } + 0{€ 



,3/2^ 



E.ng + (7/8)E,4f) 
1800 



e{VT^)\ (<0(1)) 



(3-11) 



(3-12) 



Here 'b' ('f') indicates a boson (fermion), Cg^yp|-, is as in Eq. ()2-28|) . with m = 0, 

and iYpD = l/(e'^P + 1). 

Equation (|nT|) is valid for |y| < 0(1). We note that, from Eq. the 
quantities y (oc e) and e are of the same sign. Then, from Eq. I|3-lip . we see that the 
0{^/e) term is negative and invariant under e — > — e, so that A'^be/fd ~ ^h/i is the 
same for e = |e| as for e = — |e| . The difference between the cases e > and e < 
arises at 0(e). 

Low-temperature limit: /3 — > oo 

In the limit (3p » 1, we have A^be/fd ~ e~^P, and then 



Nh/iip) 



1 - /3p + (1 + 6y)e/3p + ^(2y - l)e(/3p)2 + 0(e3/2) 



-/3p 
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The series in the square brackets on the R.H.S. here seems to be nonconvergent for 

l/V~e<Pp or l/^</3p. (3-13) 
However, due to the factor e~^^, the region (|3-13j) is unimportant. We see that 

[A^b/f(p)]e>0 < [A^b/f(p)]e = -|e|- 

High-temperature hmit: /? ^ 

In the hmiting case Pp ~ 0, we have A'be ~ ^/{Pp) and A^fd — 1/2 — Pp/A. 



Then, we obtain 



1 



Nb{p) = i- l-\/2^ + 8ye + 0(e3/2) 
f3p I 



Nfip) 



1 1 

2 ~ 4 



(3p 



1 + ^/2fe - {1 + 6y)e + 0{e^/^) 



(3-14) 



In this case, we see that [iVb/f(p)]e>o > [N^/i{p)]e=-\e\- 
3.3. Initial correlations (2n-point propagators) Cn 

Here, we compute (HlLi -^(IP«o|))c {\Pio\ = Pi), appearing in Eq. (|2-22j) . in the 
content of complex scalar field theory. A procedure similar to that above leading to 
Eq. (j3-ll|l in the present case yields [cf. Eq. H3-7p ] 

^« ( nr=i fe''^' - l) ' e2^BEV/(3/3^)' 



i=l 



2yeY,PP^{'^ + NBE{Pi)) 



i=l 



e^Ul + 6y)Ppi{l + NBE{pi)) 

i=i ^ 
2y-l 



'-{PpiYil + NBE{Pi)){l + 2NBE{Pi)) 



+ 5^ (2y - l)P'p^PJil + NBE{p^ml + Nbe{pj)) + 0(e'/') 

(3-15) 

We are now in a position to compute the connected contribution, (ITiLi -^(^'«))c- 
For the sake of generality, from this point, we consider a system of bosons and/or 
fermions. We start with n = 2, in which case we have 

{N{pi)N{p2))c = {N{p,)N{p2)) - iVb/f(pi)iVb/f(P2) . (3-16) 
Then, using Eqs. H3-ll() and (|3-15|) . we obtain 

{N{pi)N{p2))c = -el3^PlP2NBE/FD{Pl)NBE/FD{P2) (l ± iVBE/FD(Pl)) 

X (l±iVBE/FDfe)) +0(e3/2), (3.17) 
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which is of 0(e). For n = 3, we have 

= {N{pi)Nip2)Nip3)) - A^b/f(Pi)iVb/f(P2)iVb/f(P3) 

- (iV(pi)iV(p2))ciVb/f(P3) - {N{P2)N{p3))cN^/f{pi) 

- (iV(pi)iV(p3))ciVb/f(P2) 

= 0{e^/^) . (3-18) 
In a similar manner, for 4 < n, we obtain 

n n n 

(11 ^(^'^))- = (n^(^^)) - n^b/f(p*) 

1=1 1=1 i=l 

n 

- {N{pl)N{p2))cllNb/^{p^)+0{e^/^) 

perm. i=3 



Thus, we conchide that, to 0(e), C„ = for 3 < n. C2 is given by Eq. with 
n = 2 and Eq. (nTTTI) . which is of 0(e). 

Finally, we point out that in the case of Boltzmann-Gibbs statistics, we have 
(nr=i N{pi)) = m=l{N{p^)), and hence C„ = 0. 

§4. Physical implications 

4.1. Hard thermal loops in hot QCD 

According to QCD at high temperature (hot QCD), among the formally higher- 
order amplitudes are those that are of the same order of magnitude as their lowest- 
order counterparts.^^ This is the case for classes of amplitudes whose external mo- 
menta are all soft, i.e., = 0{gT) with g the QCD coupling constant. The relevant 
diagrams are the one-loop diagrams with hard loop momenta [Qioop = 0(T)], and 
for this reason, these are called hard thermal loops (HTL). Each HTL amplitude is 
proportional to the characteristic mass, the gluon thermal mass (mg), or the quark 
thermal mass (m^): 



m 



-g^T\l + Nf/6) = ^{l + Nj/6)Ibg , 



2 

3 2- ■ vr 



ml = ]:g^T^ = ^lBG, (4-1) 



where A'^^ is the number of active (massless) quark flavors. In Eq. H4-l() . we have 
introduced Xbg (where 'BG' means 'Boltzmann-Gibbs'), which is defined by 



Ibg ^ / dppNBEip) = 2 / dppNpBip) = vr^rVe . 
Jo Jo 



(4-2) 
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For Tsallis statistics, the above Ibg is modified. From Eq. (|3-11|) . we can 
compute the Tsalhs statistics counterpart of Eq. H4-2|) . Xt: 

/•oo poo 

Tt = / dppN^ip) = 2 / dppNfip) 
Jo Jo 

= 1 - 2^2^ + e(182/ - 1) + 0(e3/2)l . (4-3) 



Thus, we obtain for the gluon and quark thermal masses in Tsallis statistics, 

1 - 2v^ + e(18y - 1) + 0{e^/^) , (4-4) 



m I 

911 



\ 9liJ BG 

where ("^3/^)3(5 is as in Eq. (|4-1|) and y is as in Eq. (|3-12|) . with Udf = 16 and 
Udf = 12Nf. We find that, to 0{^/e), ^ < {ml/^)^^. 

4.2. Manifestation of the initial correlations 

Here we give some examples of quantities in which C„ participates. We consider 
a complex scalar field system. (Generalization to other field systems is straightfor- 
ward.) First, we note that each of the average values {(j)c4>c) (~ i^ccji {4'c4'c<Pc4'c) 
(~ iA^cccj, ■ ■ ■ that we know represents information that we possess concerning the 
statistical properties of the system. 

The linear response theory is formulated in terms of the so-called linear response 
function of the n-point correlation with the external field that couples to the quantity 
under consideration. The linear response functions are given by the appropriate 
Green functions, G, in the CTP formalism. As an example, we take the field 
(j)^ {(j)) itself as such a quantity. Then, the linear response function of the (n — 1)- 
point correlation (n = 2,4, •••) is given by Gc-c;c-cA iGc-A;c-c)- Similarly, the 
nonlinear response functions are also given by the appropriate Green functions in 
the CTP formalism. For example, the second-order response functions of the two- 
point correlations are given by Gcc;AA, GcA]cA and GaA;cc- 

The initial correlations, given in Eq. (|2-14p . contribute to the above-mentioned 
response functions as well as to the average correlations Gc...;c - c; and hence they 
can be "measured" in principle. For example, the diagrams that include one vertex 
and one four-point propagator C2 contribute to the four-point Green functions Gcc-m 
Gcc;cA and GcA\cc- Explicit computation of the contribution to, e.g., Gcc-cA, using 
the form for C2 obtained in §3.3 yields 

Gcc;cAixi,x2;yi,y2) 

= \ j d'^z[{2C2{z,xi\z,yi)SR{z - y2)SR{x2- z) + {xi ^ X2)} 

+C2{xi,X2] z, z)Sr{z - y2)SB{yi - z)] 



{'■ 
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+ e 



-i[Piyi~y2)+Piixi-yi)+P2ix2-yi)] 



SRiP)SR{P-Pl-P2)] 



+ {xi ^ X2} 



(4-5) 



Note that Gcc;cA{xi,X2;yi,y2) can also be rewritten in the form^^^ 
Gcc;cAixi,x2;yi,y2) 

= ^ X] [^(2^i'^2,yi,y2)( |{0(xi),(?:)(x2)} ,(/>^(yi)| ,0^(y2 

perms. 

+ 0{xi,x2,y2,yi){\^ {(I){xi),(l){x2)} ,(l)\y2) 

+0{xi,y2,X2,yi){^^ (p{xi),(p^{y2) , 0(x2)| , (/>^(yi)|) 



(4-6) 



where 'perms.' indicates that ah permutations among {xi,X2,yi) are taken, and 
9{x, y, u, v) = 9{xo—yo)6{yo—uo)6{uo—vo). From this formula, we see that Gcc-cA = 
for y2o > icio, X20, yiO) as should be the case. As must be the case, Gcc;cA in Eq. 
(|4-5|) satisfies this property. 

4.3. Remarks on higher-order corrections 

Higher-order corrections to the CTP propagators come from the following two 
sources. 

1. Corrections due to the interactions, which, as in the case of Boltzmann-Gibbs 
statistics, are obtained by computing the vertex-insertion diagrams. 

2. Corrections arising from the corrections to the pressure V computed within the 
Boltzmann-Gibbs statistics [cf. Eqs. (f^^ and (fT^ ]. For hot QCD, the 
corrections to V are known^^) up to 0[g^ ln{l/g)]. 
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